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Abstract. We study the equation ut + uu x + u — K * u = in the case of 
an arbitrary K > 0, which is a generalization of a model for radiating gas, 
in which K(y) = ^e~ l y L Using a monotone iteration scheme argument we 
establish the existence of traveling waves, which gives a solution to an open 
question raised by Denis Serre [1], 



1. Introduction 



The system of equations 



(1.1) 



{ 



u t + uu x + q x = 0, 
-qxx + q + u x = 0, 



for x G R and t > was derived by K. Hamer in 1971 [2, eqs. (4.16), (4.19)] (see 
also [3]) as a third-order approximation of the full system describing the motion 
of radiating gas in thermo-nonequilibrium (the second-order approximation there 
gives the viscous Burgers equation ut+uu x = u xx , and the first-order approximation 
gives the inviscid Burgers equation Ut + uu x = 0). The system (1.1) is formally 
equivalent to 



where K * u(x,t) = J R K(x — y)u(y,t)dy, K is even and J R K = 1. Nonlinear 
equations with a nonlocal operator being considered instead of the Laplacian have 
a long history and have recently been gaining increasing popularity, due to superior 
physical qualities of such models and at the same time often ingenious mathemat- 
ical solutions. In fact, when in 1893 J.D. van der Waals derived his variational 
thermodynamic theory of binary phase transitions, he introduced a mean field K, 
formally expanded K * w — w into power series: J R K(y)(w(z + y) — w(z))dy = 
dw"(z) + c 2 w lv + where ci = | J R y 2 K(y)dy, c 2 = ^ J R y 4 K(y)dy, etc., and 
used just the first term in this expansion [4]. He also investigated his nonlocal 
model with the Green's function K (y) = -ie '^l [4]. 

The Hamer equation (1.2) was studied in [2, 5, 6, 7, 8] and the more general 
nonlocal Burgers equation (1.3) in [9, 10, 1]. In particular, it was discovered that the 
solution of the Cauchy problem for (1.3) with smooth initial data is either smooth 
for all times or becomes discontinuous in finite time. Here we focus on traveling 



(1.2) 




and can be considered a special case of the integrodifferential equation 
(1.3) u t + uu x + u~K*u = Q, 
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wave solutions U(x — st) of (1.3). To account for possible shock discontinuities, a 
traveling wave is defined as a weak solution, i.e., a function U which satisfies 



for an arbitrary <fr £ C§° . In [6] the authors obtained the following result (see also 
a partial version in the earlier [5]). 



Proposition 1.1. [6] For every pair u_ > u + there exists a solution U(x — st) 
of (1-4) with K{y) — ^e^^ , U(—oo) = U-, ?7(+oo) = u + and s = \{u- +u + ). 



u_ - u+ < V2, U is C 1 . (c) When u- - u+ < U is C k . 

Subsequently, Denis Serre raised the question whether a similar result can be 
established for (1.3) [1, Open Problem 1, p. 33]. In Section 2 we answer it in the 
affirmative, though we fall short of providing an exact regularity criterion of the 
profiles as in Proposition 1.1. 



Theorem 2.1. Let K >0,K& W 1 ' 1 (R), K nonincreasing on (0,oo), f R y 2 K(y)dy 
< oo. For every pair u_ > u+, there exists a solution U(x — st) of (1.4) such that 
U'(x) < for x 7^ 0, U{—oo) = U-, U(+co) = u + and s = ^(u_ +u + ). Moreover, 
if U- — u + > 4 J R \y\K(y)dy, then U is discontinuous at one point. 

Proof. Proposition 1.1 was proved by a phase-plane analysis [6]. Here we argue 
in an indirect way as we construct a monotone iteration scheme and pass to the 
limit. This is a standard elliptic technique, sometimes called the sub-supersolution 
method, whose use in nonlocal problems dates to [11] (see also the many follow- 
up papers on traveling waves in various nonlocal monostable models). In an early 
and unpublished version of [12] the authors used it to construct continuous and 
discontinuous stationary waves of the nonlocal bistable equation [13]. A similar 
argument appeared in [14] for the bounded interval case. 

The main difficulty is to set up the difference scheme in a way that the sequence 
of monotone solutions is monotone. The whole sequence then has a limit so that 
one can pass to the limit on both sides of the scheme. This procedure would not 
work if only a subsequence was convergent. To ensure the limit is not trivial, one 
also needs to find an appropriate subsolution. 

It suffices to show that for any constant u c > 0, 

(2.1) uu' = K * u - u 

admits a monotone and odd solution u such that u(— oo) = u c and m(+oo) = — u c . 
If we let U = u + s, where u c — \{u- — u + ) and s = \{u- + u + ) (the Rankine- 
Hugoniot condition), then, since J R K = 1, U satisfies 



(1.4) 





2. Existence of traveling waves 



U'(U-s) =K*U-U, 



U(—oo) = it_ and U(+oo) = u + . 
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We will see that it suffices to study (2.1) on (— oo,0). We find it convenient first 
to find a supersolution and a subsolution for (2.1). Let 

!u c , x < 0, 
0, x = 0, 
— u c , x > 0. 

Then it is easily seen that 

MO + UqUq > K * Uq 

on (—oo,0), thus u a is a supersolution there. 
Let 

s(x) = -tan -1 (err). 

We investigate the function 

- f K(x-y) (tan- 1 (ey) -tan- 1 (ex)) dy 
g(x,e) = k • 

tan (ex)- ^—^ 

7T 1 + e 2 x 2 

Using PHopital's rule, we get lim a; ^_ 00 g(x, e) = 0, lim £ ^ g( x , e ) = for any x < 0, 
and 

Um <A:, «) - -li^l±gi±Zlj±l^ . « / £[M 
^o yv 7 2w c ^_ 2xe 3 2u c J R l + e 2 y 2 y 

1 + e 2 .T 2 1 + 6 2 ir 2 ^ (1 + e 2 x 2 ) 2 

Since J fl y 2 K(y)dy < oo, we see that for small enough e > 0, e) < 1 for x < 0, 
therefore s satisfies 

s + ss 1 < K * s, 
thus s is a subsolution on (— oo, 0). 

On (— oo, 0) define the sequence {u n } n >o by 

(2.2) u n+ i + u n u' n+1 = K *u n 

and u n (— oo) = u c . This iteration has the same properties on (— oo,0) and (0,oo), 
by odd extension and the evenness of K. In the inductive step we show that the 
statement 3 continuous u n such that u n > 0, u n (— oo) = u c , u' n < 0, u n — < 
0, uo > u n > s implies the statement 3 continuous u n+1 such that u n+1 > 0, 
u n+ i(-oo) = u c , u' n+1 < 0, u n+ i - u n < 0, u > u n+ i > s. 

The solution to the linear equation (2.2) with the initial condition tt„ + i(0) = 
u n+i > an( i ^ < is given by the formula 

, n f ds r n f f° K * U„ (s) 

(2.3) u n+1 (x) = e'- ^) [u° n+1 - / e" ^> - 

However, u n+ i in (2.3) is such that u n+ i(— oo) = u c only if we choose 

\ n /"° ro dr K*U n (.<i) 

which follows from l'Hopital's rule: 



,o +1 _ [°e-f^ K * U ^ d s 

J X 



lim ^ ^5-^ = lim X * u„(z) = u c 

x^-cc e J x ^J^J x^-oo 
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_ f° dr s 

Note that the integral in (2.4) is convergent, since u n < u c implies e Js < e u <= . 
Let v n+1 = u n+x - u n , v n = u n - u„_i. From (2.2) we get 

Vn+l + U n u' n+1 - U n -iu' n = K * V n . 

After adding and subtracting u n u' n on the left side we get 

(2.5) v n+1 +u n v' n+1 +u' n v n = K *v n . 
Since v n is odd, K * v n can be represented on (— oo, 0) as 

(2.6) K*v n (x) = f [K(x-y)-K(x + y)]v n {y)dy. 

J —oo 

Notice that K(x — y) — K(x + y) > for x < and y < because K is even and 
non-decreasing for x < 0, thus K * v n > 0. If v n+ \ > somewhere, then since 
v n +i{— oo) = 0, 3^0 such that v n +i(xo) > and v' n+1 (xQ) > 0. Since v n < 0, 
u ' n < 0, the left side of (2.5) is positive, a contradiction. Now differentiate (2.2): 

(2.7) u' n+1 + u' n u' n+1 +u n u' n+1 = K'*u n = -K{x)[u n (0-)-u n ({)+)] + K *u' n < 0. 

If u' n+1 > somewhere, then since u n+ i < u n , 3xq such that u n+1 (xo) > and 
Un + i(xo) > 0. The left side of (2.7) at x is then positive since u' n (x ) > —1 from 
(2.2), a contradiction. Using u' n+1 < Owe get u n+ \ > u n+ \ +u n u' n+1 = K *u n > 
for x < 0. 

Assuming u n < u we get 

u n +i + u u n+1 < u n+ i + u n u' n+l = K *u n < K *u =u 1 + u u[, 

thus 

u n+ i - ui + u (u' n+1 - u[) < 0, 

and since u n+ i(— oo) = Uo(— oo), u n+ \ < u by a similar argument as above. 
Assuming u n > s, in the same way we get 

u n +i - ui + s(u' n+1 - u[) > 

and > s. 

Thus {u n } is a sequence of monotone functions such that u n+ i(x) < u n (x) for 
all n > and x < 0. Let u = linin^ooUn. Since uo > u n > s for all n > 0, 
u(— oo) = u c . In (2.2) add and subtract u n +iu' n+1 , multiply it by cf) e C°°, with ^ 
compactly supported in (— oo,0), and integrate over (— oo,0): 

(2.8) / [--u n+1 (f> + u n+ i4> - u' n+l {u n+ i - u n )4> - (K * u„)4>] = 0. 

J — oo 4 

From uq > u n > s and (2.2), \u' n \ is uniformly bounded in n, thus we can use the 
Lebesgue Dominated Convergence Theorem in (2.8) and pass to the limit as n — > oo 
to get 



r 1 

/ [-^u 2 (f>' + (u- K *u)(f>] =0. 



Since \u 2 G VF 1 '°°(— oo, 0), ^u 2 is diffcrentiable a.e. and {\u 2 )' = K * u — u a.e. 
on (— oo, 0). \v? G T / F 1,00 (— oo, 0) also implies that and therefore u is Lipschitz 
continuous, thus u G C 1 (— oo,0) and u satihes (2.1) on (— oo,0). Since u is odd, it 
also satisfies (2.1) on R\ and (1.4). U'(x) < for x ^ is shown in a similar way 
as the argument below (2.7). 
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To find a condition for which u is discontinuous, we argue by contradiction. 
Assume that u is continuous at 0. Note that v! is integrablc on R, since for x ^ 
we have u' — — 1 > —1. Integrating (2.1) over (— oo,0) and using Fubini's 
Theorem and then the LDC Theorem we get 

— -u c 2 = f uu' = f (K*u — u)=f f K{y) I yu' (x + yt)dtdydx 

^ J-oo J-oo J-ooJR Jo 

= / yK(y) / lim [u(yt) - u{-m + yt)]dtdy = / yK(y) \ u(yt)dtdy. 
Jr Jo J R Jo 

The right side can be estimated to yield 

i« c 2 < u c J \y\K(y)dy. 

Thus if u c 2 > 2 J R \y\K(y)dy we obtain a contradiction. This gives the criterion for 
a solution to be discontinuous as it_ — u + > 4 J R \y\K(y)dy, since u c = \{ u - — u +)- 
For K{y) = | e - fc ^, we get U- - u+ > |. □ 
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